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ABSTRACT 

Progesterone receptors (PR) are transcription 
factors relevant to breast cancer biology. Herein, 
we describe an N-terminal common docking (CD) 
domain in PR-B, a motif first described in mitogen- 
activated protein kinases. Binding studies revealed 
PR-B interacts with dual-specificity phosphatase 6 
(DUSP6) via the CD domain. Mutation of the PR-B 
CD domain (mCD) attenuated cell cycle progres- 
sion and expression of PR-B target genes 
(including STAT5A and Wnt1); mCD PR-B failed 
to undergo phosphorylation on Ser81, a ck2- 
dependent site required for expression of these 
genes. PR-B Ser81 phosphorylation was depend- 
ent on binding with DUSP6 and required for re- 
cruitment of a transcriptional complex consisting 
of PR-B, DUSP6 and ck2 to an enhancer region 
upstream of the Wnt1 promoter. STAT5 was 
present at this site in the absence or presence 
of progestin. Furthermore, phospho-Ser81 PR-B 
was recruited to the STAT5A gene upon progestin 
treatment, suggestive of a feed-forward mechan- 
ism. Inhibition of JAK/STAT-signaling blocked 
progestin-induced STAT5A and Wnt1 expression. 
Our studies show that DUSP6 serves as a 
scaffold for ck2-dependent PR-B Ser81 phos- 
phorylation and subsequent PR-B-specific gene 
selection in coordination with STAT5. 
Coregulation of select target genes by PR-B and 
STAT5 is likely a global mechanism required for 
growth promoting programs relevant to mammary 
stem cell biology and cancer. 



INTRODUCTION 

Progesterone is an ovarian steroid hormone essential for 
breast development and implicated in breast cancer pro- 
gression. Progesterone receptors (PR) exist primarily as 
two coexpressed isoforms, PR-A and PR-B (1,2), 
encoded by the same gene downstream of distinct pro- 
moters (3). PR-B, the full-length receptor, contains 164 
amino acids at the N-terminus, not present in PR-A, 
termed the B-upstream segment (BUS) (4). Both receptors 
contain the same DNA-binding domain (DBD), a hinge 
region (H) and two activator function (AF) domains; PR- 
B contains a third AF domain in the BUS (5). Unliganded 
PR rapidly shuttles between the cytoplasm and the 
nucleus. After ligand binding, however, PR undergoes 
dimerization and is retained in the nucleus. Nuclear PR, 
together with coactivators and corepressors, activates or 
represses transcription of PR target genes, either directly 
through DNA binding to progesterone response elements 
(PREs) or indirectly through tethering interactions with 
other transcription factors (API, SP1, STATs) (6-9). 

PR-mediated regulation of gene expression is controlled 
by many posttranslational modifications to the receptor, 
primarily on N-terminal serine (phosphorylation) and 
lysine (ubiquitination, sumoylation and acetylation) 
residues (10-15). These modifications significantly alter 
receptor stability, localization, transcriptional activity 
and target gene selectivity. PR is phosphorylated on 
serines (Sers) 294, 345 and 400 by mitogen-activated 
protein kinase (MAPK) and cyclin-dependent kinase 2 
(cdk2). PR-B is also phosphorylated on Ser81 by ck2 [for- 
merly casein kinase II; (16-19)], a ubiquitously expressed, 
constitutively active kinase that is overexpressed in every 
cancer examined thus far, including breast cancer (16,20). 
ck2-dependent PR-B phosphorylation of Ser8 1 regulates a 
specific subset of PR-B target genes involved in breast 
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cancer cell growth and pro-survival, including BIRC3, 
HSD11(32 and HbEGF (19). In addition, ck2 is recruited 
along with Ser81-phosphorylated PR-B to enhancer sites 
of a subset of progesterone-responsive target genes (19). 
Notably, these studies have shown that ck2-dependent 
phosphorylation of PR-B Ser81 is unique to PR-B and 
thereby a primary determinant of PR isoform-specific 
activity. However, the molecular interactions necessary 
to support PR-B Ser81 phosphorylation have yet to be 
understood. 

Posttranslational modifications to PR regulate the recep- 
tor's interactions with other proteins (21). PR protein inter- 
action domains include the estrogen receptor (ER) 
interaction domains [ERIDs (22)] and a poly-proline-rich 
(p-Pro) domain that is required for interaction between PR 
and the SH3-domain of c-Src (23). PR interacts with many 
other proteins via unknown domains (i.e. MEK1, FGFR2, 
STAT5) (21,24-26). In a recent in silico analysis of the PR 
amino acid sequence aimed at identifying protein inter- 
action domains, we identified a putative common 
docking (CD) domain in the N-terminal BUS region of 
full-length PR-B (21), a region that is not present in other 
PR isoforms. CD domains are commonly found in 
members of the MAPK family, where they mediate 
interactions between MAPKs and their upstream activa- 
tors (MEKs), negative regulators (MAPK phosphatases 
[MKPs] or dual-specificity phosphatases [DUSPs]) and 
downstream targets (27,28). CD domains are characterized 
by clusters of negatively charged amino acids (aspartic or 
glutamic acid) that form electrostatic interactions with a 
positively charged 'D domain' in their respective binding 
partner [MKKs, DUSPs, MAPK substrates; reviewed 
in (29)]. The PR-B CD domain is an identical match to 
the CD domain of the MAPK family member, Erk2. This 
identical match suggests that PR-B interacts with the 
same D domain-containing proteins as Erk2. However, 
the function of this domain, unique to the PR-B isoform, 
has not yet been determined. 

We predict that both PR-B Ser81 phosphorylation and 
PR-B CD domain interactions may be involved in breast 
cancer progression. Indeed, PR has been implicated in 
breast cancer progression in recent clinical studies of 
hormone-replacement therapy. These studies found that 
women taking estrogen and progesterone had more and 
larger breast tumors than women taking estrogen 
alone (30,31). In addition, we recently identified a 
phosphorylated PR-B gene signature associated with 
decreased survival in women with luminal (ER+ and/or 
PR+) breast cancer whose disease stopped responding to 
anti-ER therapy with tamoxifen (13). Mitogenic protein 
kinases — such as MAPK, c-Src, cdk2 and ck2 — are 
commonly upregulated in cancer (17,32-34) and represent 
likely pathway components in PR-mediated gene expres- 
sion in breast cancer. Understanding how PR interacts with 
these protein kinases and their regulatory protein partners 
is critical to understanding breast tumor etiology and de- 
veloping better treatments. Herein, we sought to identify 
proteins that interact with PR-B via a novel N-terminal CD 
domain and how protein-protein interactions through this 
domain alter PR-B phosphorylation and transcription 
factor function in breast cancer models. 



MATERIALS AND METHODS 

Cell lines and constructs 

The estrogen-independent ER/PR-positive T47Dco 
(T47D) variant cell line has been previously described 
(35), and is the parent cell line from which all T47D- 
variants used herein were created. T47D-Y (PR negative), 
T47D-YB (stably expressing wt PR-B) and T47D-YA 
(stably expressing wt PR-A) cells were characterized by 
Sartorius et al. (36). T47D-S79/81A PR-B cells have been 
previously described (19). T47D-mCD PR-B cells were 
created by stable expression of pSG5-mCD PR-B and 
pSV-neo in T47D-Y cells using FuGene-HD (Roche). 
Individual colonies were selected in 500ug/ml G418 and 
maintained in 200ug/ml G418 after initial selection in 
cMEM (as described in the Supplementary Materials and 
Methods). All cells were maintained as previously 
described or described in the Supplementary Materials 
and Methods (19). Creation of the pSG5-mCD PR-B, 
CD-PR-A and mCD-PR-A plasmids was performed as 
described in the Supplemental Materials and Methods. 
The DUSP6 construct [pcDNA3.1(-)Myc-His] was a gift 
from Stefanie Dimmeler, University of Frankfurt. 

Gene expression profiling 

T47D cells stably expressing pSG5 empty vector, wt PR-B 
or mCD PR-B were serum starved in modified iMEM 
(Gibco) for 1 day and treated with R5020 (10 nM) or 
vehicle control for 6h. Total RNA was extracted using a 
RNeasy kit (QIAgen) with on-column DNase I treatment 
(QIAgen). Triplicate RNA samples were labeled and 
hybridized to the Illumina HT-12v4 bead chip platform 
according to the manufacturer's protocols. Chip scanning 
within Genome Studio software produced raw expression 
values that were analyzed within R software using the 
Bioconductor (37) package called lumi, where raw 
intensities were log2 transformed and quantile normalized. 
Differentially expressed genes were analyzed using the 
limma package, where empirical Bayes was used to 
better estimate the variance of the genes. Reported gene 
expression data contain log2 normalized intensities and 
biological comparisons (e.g. R5020/vehicle) contain log2 
fold change with the Benjamini and Hochberg (BH)- 
adjusted P-value (38). Heat maps were generated 
through unsupervised hierarchical clustering of probes 
by the heatmap.2 function in the gplots R software 
package. Clustering was performed using Euclidean 
distance and complete linkage. Rows were scaled to 
have mean zero and standard deviation (SD) equal to 
one. Gene expression data are available in the NCBI 
Gene Expression Omnibus database (accession number: 
GSE46850. 

Pathway and gene set analysis 

Ingenuity Pathway Analysis (IPA) software (Ingenuity® 
Systems, www.ingenuity.com) was used to compare biolo- 
gical functions or network pathways in cells expressing 
wt or mCD PR-B after progestin treatment (wt+R5020/ 
wt - R5020 or mCD + R5020/mCD - R5020). Default IPA 
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settings for Core analyses were used to evaluate 
upregulated genes (fold change >1.5, P<0.05). IP A 
Comparison analyses were used to reveal whether cells ex- 
pressing wt or mCD PR-B regulated functionally distinct 
pathways. Analyses were scored based on significance (the 
BH-adjusted P-value, corrected for multiple hypothesis 
testing) and the threshold for a gene list to be significantly 
involved in a particular biological function was P < 0.05 
[i.e. -log 10 (BH-adjusted P-value) >1.30]. 

Gene set enrichment (GSEA) and Leading Edge 
analysis (39,40) was performed using the javaGSEA 
desktop software; all five gene set collections (cl-c5) 
from the Molecular Signatures Database (MSigDB) 
version 3.1 were queried. Data set files were developed 
based on normalized Illumina expression intensities from 
cells that constitutively express PR. Specifically, the log2 
fold change values were compared for two phenotypes in 
our ligand-dependent analysis: (wt +R5020/-R5020) 
versus (mCD +R5020/-R5020). GSEA was executed 
using the default settings, except the permutation type 
was set to Gene_set with 1000 permutations, and the 
metric for ranking genes was set to Diff_of_Classes 
because normalized expression data were log2 trans- 
formed. Each MSigDB collection was analyzed individu- 
ally in multiple GSEA runs. 

Immunoblotting 

Immunoblotting was performed as previously described 
(19) and in the Supplementary Materials and Methods. 

Coimmunoprecipitation experiments 

For coimmunoprecipitation (CoIP) experiments, cell 
lysates were collected in RIPA (supplemented as described 
in the Supplementary Materials and Methods) and 
incubated on ice for 60min. Cell lysates containing equiva- 
lent protein concentrations (1000 ug) were incubated over- 
night at 4°C with 2 ug appropriate antibody or control IgG. 
Protein G agarose (Roche Diagnostics) was added for the 
final 1 h of incubation time. Immune complexes were 
washed three times with supplemented RIPA buffer, 
resuspended in Laemmli sample buffer containing 
dithiothreitol and P-mercaptoethanol, boiled for 5 min 
and subjected to western blotting analysis. 

Transient transfections 

Transient transfections were performed as previously 
described (19) and in the Supplementary Materials and 
Methods. 

H 2 0 2 assays 

After T47D-wt PR-B cells were starved overnight in 
serum-fee iMEM, they were pretreated with 1 mM H 2 0 2 
for 20 min, followed by 30 min with 10 nM R5020. Protein 
lysates were isolated and analyzed as described above. 

Luciferase transcription assays 

Luciferase assays were performed as previously described 
(41) using the Dual Luciferase Reporter Assay (Promega). 
Relative luciferase units were normalized to Renilla ± SD. 



siRNA 

ON-TARGETplus SMARTpool designed to target human 
DUSP6 (DUSP6) and nonsilencing siRNA controls were 
purchased from Dharmacon. For siRNA experiments, 
T47D-YB cells were plated, and 24 h later they were trans- 
fected with 50 nM nonsilencing or DUSP6 siRNA. After 
72 h, cells were treated with EtOH or lOnM R5020 for 
60 min. Protein lysates were isolated and analyzed as 
described in the Supplementary Materials and Methods. 

Reagents 

Cells were treated with the following reagents (when 
applicable): R5020 (lOnM; Sigma), AG490 (50 uM; 
CalBioChem) and H 2 0 2 (ImM; from Cell Biolabs 
OxiSelect™ Intracellular ROS Assay Kit). 

Cell cycle analysis/flow cytometry 

Flow cytometry for cell cycle analysis was performed as 
previously described (42). 

Real-Time Quantitative polymerase chain reaction 

Real-time quantitative polymerase chain reaction (RT- 
qPCR) was performed as described previously (19) and 
in the Supplementary Materials and Methods. Primer 
sets used for qPCR are listed in Supplementary Table 
S2. Relevant genomic sequence information and 
enhancer positioning for Wntl genomic sequence is 
based on GR37 Release 57 (March 2010). 

ChIP assays 

ChIP was performed using the ChlP-IT Express Kit 
(Active Motif) according to the manufacturer's instruc- 
tions using sonication as the method for chromatin 
shearing. Lysates were immunoprecipitated (IP) overnight 
(18 h) with the following antibodies: PR (ThermoScientific 
#MS-298-P), STAT5 (Santa Cruz Biotechnology sc-1081 
and sc-836), ck2a (Santa Cruz sc- 12738), DUSP6 (Santa 
Cruz Biotechnology sc- 100374) or an equal amount of 
mouse or rabbit IgG. Resulting DNA was analyzed 
using qPCR as described above, and data are represented 
as a percentage of input DNA. In silico analysis using 
Matlnspector (Genomatix) identified potential PRE- 
binding sites using the following consensus sequence: 
RGNACANRNTGTNCY. Primer sets used for ChlP- 
qPCR are listed in Supplementary Table S2. 

CEAS 

Web-based CEAS analysis (http://ceas.cbi.pku.edu.cn/ 
index.html) was performed on a publicly available PR 
ChlP-Chip data set (http://cistrome.org/NR_Cistrome/ 
Cistrome.html). TRANSFAC and JASPAR motifs 
were used to determine putative transcription factor 
binding sites. 

Statistics 

Statistical significance for all experiments was deter- 
mined using an unpaired Student's /-test, unless otherwise 
specified. 



Nucleic Acids Research, 2013, Vol. 41, No. 19 8929 



RESULTS 

PR-B CD domain is required for progestin-induced 
S-phase entry 

We previously identified a putative CD domain 
(Figure 1A) (21) located in the N-terminal BUS region 
of full-length PR-B, a region that is absent from other 
PR isoforms (Figure IB). To study the importance of 
this newly identified CD domain in modulating PR-B- 
specific functions, we mutated the critical negatively 
charged amino acids (D-aspartic acid, E-glutamic acid) 
to alanines (A), creating an mCD PR-B mutant 
(Figure IB). Similar mutational strategies have been 
used to study CD domains present in Erks and other 
MAPKs (27,43). We then determined whether mCD PR- 
B was able to bind DNA and activate PRE-dependent 
transcription in luciferase reporter gene assays. PRE- 
luciferase expression levels were increased at similar 
levels in HeLa cells transiently expressing wt or mCD 
PR-B following treatment with vehicle (ethanol) or 
10 nM R5020 (synthetic progesterone) (Figure 1C). 
These data suggest that the PR-B CD domain is not 
required for intrinsic PR-B transcriptional activity. 

Progesterone or synthetic progestins (R5020) induce S- 
phase entry in breast cancer cells expressing PR-B (44^16), 
but not PR-A (47). Experimental isolation of PR isoform- 
specific actions is complicated by the fact that estradiol is 
generally required for robust PR expression in steroid 
hormone receptor-positive breast cancer models. Not 
only is estrogen itself a potent mitogen, but it tightly 
controls PR isoform expression (48). To overcome this 
barrier, we used the ER-positive T47Dco cell line, which 
expresses abundant PR-A and PR-B in the absence of 
added estrogen (35). A naturally occurring PR-negative 
variant of T47Dco cells, termed T47D-Y, was used to 
create stable cell lines constitutively expressing either wt 
PR-B (T47D-YB) or wt PR-A (T47D-YA) (49). We then 
tested the contribution of the PR-B CD domain to cell 
cycle progression by stably expressing mCD PR-B 
(T47D-mCD PR-B) in T47D-Y cells and analyzed proges- 
tin-induced S-phase entry. As predicted, there was an 
increase in progestin-induced S-phase entry in T47D-YB 
cells but not in T47D-YA cells (Figure ID). Interestingly, 
cells stably expressing mCD PR-B also failed to enter 
S-phase upon treatment with R5020; these cells resembled 
PR-A expressing cells (Figure ID). These data suggest 
that the PR-B CD domain is essential for proliferative 
signaling in breast cancer cells, as measured by proges- 
tin-induced S-phase entry. 

PR-B CD domain regulates select PR-B target genes 

Although mutation of the PR-B CD domain did not ap- 
preciably alter the absolute levels of PR-B transcriptional 
activity, posttranslational modifications can dramatically 
alter PR target gene selectivity, directing PR to specific 
enhancer or promoter regions in chromatin (50-52). This 
is an important feature of steroid hormone receptor action 
that is missed in reporter assay systems such as luciferase. 
To determine whether the PR-B CD domain functions 
in the regulation of endogenous PR target genes, we 
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Figure 1. CD domain in PR-B is required for progestin-induced 
S-phase entry. (A) Table comparing the CD domain core amino acid 
sequences in PR-B and selected MAPKs. (B) Schematic representing 
functional domains in wt PR-B. The full-length receptor has an 
N-terminal region unique to PR-B, termed the BUS. Both PR-B and 
the shorter form, PR-A, contain two activating function (AF 1 and 2) 
domains, a DNA-binding domain (DBD), hinge region (H) and a 
hormone-binding domain (HBD). PR-B contains an additional AF3, 
located in the BUS. Arrows mark the translational start sites for PR-B 
and PR-A. The CD domain is located between amino acids 68-76, 
immediately upstream of a ck2-dependent phosphorylation site, 
Ser81. Negatively charged amino acids were mutated to alanines (in 
red) to create the mCD PR-B mutant. (C) Reporter gene assays 
measuring transcriptional activity of PR constructs. HeLa cells were 
transiently transfected with plasmids expressing wt PR-B, mCD PR-B 
or vector, as well as a firefly PRE-luciferase reporter construct and 
Renilla expression control. Luciferase assays were performed as 
described in the 'Materials and Methods' section. Fold relative 
luciferase units (PRE-luciferase over Renilla luciferase controls) of 
R5020/EtOFf-treated cells are plotted. Error bars represent ± SD of 
three independent experiments. (D) PR-B CD domain required for pro- 
gestin-induced S-phase entry. T47D cells stably expressing wt PR-B 
(T47D-YB), PR-A (T47D-YA) or mCD PR-B (T47D-mCD PR-B) 
were starved for 18 h in serum-free media, followed by treatment 
with 10 nM R5020 or vehicle for 18h. Single cells were analyzed by 
flow cytometry. 
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performed global gene expression analyses using Illumina 
HT-12v4 whole genome bead arrays. Triplicate gene ex- 
pression analyses were performed on T47D-Y (PR-null), 
T47D-YB and T47D-mCD PR-B cells after 6 h of treat- 
ment with R5020 or vehicle (ethanol). Transcriptional dif- 
ferences between cells expressing wt and mCD PR-B 
are evident in the heat map of significantly upregulated 
or downregulated (> 1.5-fold) genes (Figure 2A). 
Differential regulation of numerous genes that require 
an intact CD domain for ligand-induced expression was 
validated by RT-qPCR (Figure 2B). Using IPA software, 
we compared the wt and mCD PR-B gene sets to a large 
database of genes that have been manually assigned to 
molecularly defined pathways, biological functions or 
disease states. Interestingly, genes that were specifically 
upregulated (> 1.5-fold) in cells expressing wt but not 
mCD PR-B (P < 0.05, BH adjusted) were identified as sig- 
nificantly involved in pathways regulating cell prolifer- 
ation, survival and cancer (Figure 2C; line represents 
P = 0.05 significance). Importantly, the individual genes 
validated by RT-qPCR (Figure 2B) were included in the 
CD-regulated gene set assigned to these IPA-defined 
pathways. These data suggest that the CD domain in 
PR-B is critical for PR's ligand-dependent contributions 
to cell growth and survival pathways, and verify that the 
CD domain regulates a biologically significant subset of 
PR-B target genes. 

PR-B CD domain is required for PR-B Ser81 
phosphorylation in response to ligand 

PR phosphorylation occurs on multiple sites and is a key 
determinant of receptor localization, ubiquitin-dependent 
turnover, tethering interactions and hormone responsive- 
ness at selected PR target genes [reviewed in (53)]. We 
therefore screened for differences in basal and regulated 
phosphorylation of wt and mCD PR-B using phospho- 
specific PR antibodies. Notably, Ser81, a basally 
phosphorylated site that is further upregulated by ck2 in 
response to ligand-binding (19), failed to undergo basal 
(absence of ligand) phosphorylation in HeLa cells transi- 
ently expressing mCD PR-B (Figure 3A) or T47D cells 
stably expressing mCD PR-B (Figure 3B). Similarly, 
ligand-induced PR-B Ser81 phosphorylation was greatly 
diminished in cells expressing mCD PR-B relative to wt 
PR-B (Figure 3 A and B). A time course of PR-B Ser81 
phosphorylation in response to R5020 treatment verified 
that mCD PR-B is persistently weakly phosphorylated 
relative to wt PR-B (Figure 3C). Additionally, we 
analyzed PR phosphorylation on proline-directed sites 
(Sers 294, 345 and 400) in HeLa cells transiently trans- 
fected with either wt or mCD PR-B and then treated 
with R5020. Despite relatively equal levels of total wt or 
mCD PR-B expression, PR phosphorylation occurred 
with faster kinetics in cells expressing mCD PR-B 
relative to cells expressing wt PR-B. After 60min, 
however, similar levels of phosphorylation were achieved 
in both groups (Supplementary Figure SI). These data 
suggest that mutation of PR-Bs CD domain dramatically 
alters PR phosphorylation in response to ligand. 
Namely, Ser81 fails to be persistently phosphorylated, 
while a number of proline-directed sites (Sers 294, 345 
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Figure 2. Gene expression analysis in mCD PR-B cells. (A) Heat map 
highlighting transcriptional differences between cells stably expressing 
wt PR-B, mCD PR-B or PR-null. The indicated T47D cells were 
treated with R5020 or vehicle. Genes differentially expressed >1.5 
fold (BH-adjusted /><0.01) are displayed for each treatment group. 
The experiment was performed in triplicate. (B) Validation of CD- 
regulated genes. T47D-Y cells stably expressing either wt or mCD 
PR-B were starved for 18 h in serum-free media, followed by treatment 
with 10 nM R5020 or EtOH for 6h. mRNA levels for selected genes 
were analyzed by qPCR. Error bars represent ± SD. (C) IPA 
comparing > 1.5-fold upregulated genes from wt PR-B and mCD PR- 
B expressing cells. Line represents P = 0.05 significance. 
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and 400) appear to exhibit transient or temporary hyper- 
phosphorylation that is not persistently maintained. 

PR-B CD domain interacts with DUSP6 

The fact that mCD PR-B lacks Ser81 phosphorylation 
suggests that the CD domain may facilitate a specific 
interaction(s) between PR-B and one or more factors 
that are required for this phosphorylation event. An inter- 
action between ck2 and DUSP6 [also called MAPK 



phosphatase 3 (MKP3)] has previously been reported 
(54). To test whether PR-B also interacts with DUSP6, 
COS cells were transiently cotransfected with constructs 
encoding wt or mCD PR-B and DUSP6 (myc-tagged) or 
vector-only controls. In cells transfected with wt PR-B 
and DUSP6, myc-tagged DUSP6 clearly coimmuno- 
precipitated with wt PR-B in a largely ligand-independent 
manner (Figure 4A; lanes 3^1). In contrast, myc-tagged 
DUSP6 weakly copurified with mCD PR-B under the 
same experimental conditions (Figure 4A; lanes 7-8), 
indicating that the CD domain is mediating an interaction 
between PR-B and DUSP6. 

To determine the specificity of PR-Bs interaction with 
DUSP6, we engineered two additional PR-specific 
mutants (Figure 4B) in which a wt or mutant PR CD 
domain was fused to the N-terminus of PR-A (termed 
CD-PR-A and mCD-PR-A, respectively). COS cells were 
transiently transfected with control (wt PR-B) and mutant 
PR constructs, as well as myc-tagged DUSP6. 
Reproducibly, wt PR-B and DUSP6 exhibited robust 
interaction as measured by CoIP, and mCD PR-B again 
displayed greatly reduced interaction with DUSP6 
(Figure 4C; lanes 2 and 4). Wt PR-A coimmuno- 
precipitated with low levels of myc-tagged DUSP6, 
similar to levels observed for mCD PR-B (Figure 4C; 
lanes 5 and 4). Transient expression of CD-PR-A and 
mCD-PR-A fusion proteins remained poor relative to wt 
PR-A (Figure 4C; lanes 6-7 of Input lysates). However, wt 
PR-A and both PR-A fusion proteins were visible in 
western blots of immunoprecipitates (Figure 4C; lanes 
5-7 of IP:PR blot). Despite relatively poor expression, 
the CD-PR-A fusion receptor coimmunoprecipitated 
with myc-tagged DUSP6 at greater levels than wt PR-A 
(Figure 4C; lanes 6 and 5). However, the mCD PR-A 
fusion receptor completely reversed this effect 
(Figure 4C; lane 7), returning DUSP6 CoIP levels to ap- 
proximately those seen with wt PR-A alone (no CD 
domain) and mCD PR-B. An additional region(s) 
common to both PR-B and PR-A (i.e. outside the BUS) 
is likely capable of weak interaction with DUSP6 (lanes 5 
and 7). These data indicate that the PR-B CD domain 
is primarily responsible for mediating the interaction 
between wt PR-B and DUSP6. 

DUSP6 is required for ck2-dependent PR-B Ser81 
phosphorylation 

We next sought to determine how PR-B's interaction with 
DUSP6 is related to PR-B Ser81 phosphorylation. We 
previously identified PR-B Ser81 as a ck2-dependent site 
regulated in response to treatment of breast cancer cells 
with progestin, and during the S-phase of the cell cycle in 
the absence of progestin (19). If DUSP6 primarily func- 
tions to recruit ck2 for PR-B Ser81 phosphorylation, then 
loss of DUSP6 should block this phosphorylation event. 
To test this hypothesis, a DUSP6-specific siRNA was used 
to knock down DUSP6 protein expression in breast cancer 
cells before analysis of progestin-induced PR-B Ser81 
phosphorylation. Although DUSP6 knockdown efficiency 
remained weak (~50% as determined by densitometry), 
T47D-YB cells transfected with DUSP6 siRNA 
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Figure 4. PR-B interacts with DUSP6 through the CD domain. (A) 
CD-domain dependent interaction between PR-B and DUSP6. COS 
cells were transiently transfected with wt PR-B, mCD PR-B or vector 
alone, as well as myc-tagged DUSP6 (or vector). Twenty-four hours 
after transfection, cells were starved for 1 8 h in serum-free media and 
then treated with 10 nM R5020 or EtOH for 60min. Lysates were IP'd 
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(25 ug not subjected to immunoprecipitation) were analyzed via western 
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antibodies. Myc-tagged DUSP6 resolves as a doublet by western 
blotting due to alternate translational start sites. (B) Schematic repre- 
senting wt and mutant PR-B constructs. The CD domain in wt PR-B 
(top) was mutated to create mCD PR-B (as described above). PR-A 
lacks a BUS region where the CD domain is located. A synthetic wt or 
mutant CD domain was fused to the N-terminus of wt PR-A to create 
CD-PR-A or mCD-PR-A, respectively. (C) CD domain-mediated 
interaction between PR-B and DUSP6. COS cells were transiently 
transfected with the indicated PR constructs or vector alone, as well 
as myc-tagged DUSP6 (or vector). CoIP and western blotting were 
performed as described in (A). 



consistently exhibited decreased PR-B Ser81 phosphoryl- 
ation relative to cells transfected with nonsilencing control 
siRNA (Figure 5A, left panel); a 50% decrease in DUSP6 
protein levels resulted in at least 75% reduction of PR-B 
Ser81 phosphorylation (Figure 5A, right panel). As a 
control for functional DUSP6 knockdown, we measured 
Erkl/2 phosphorylation under similar conditions because 
DUSP6 phosphatase activity is a negative regulator of 
Erkl/2 phosphorylation. As expected, T47D-YB cells 
transfected with DUSP6 siRNA contained elevated 
levels of Erkl/2 phosphorylation relative to controls, 
indicating effective DUSP6 knockdown (Figure 5B). To 
confirm these results using independent methods, we 
chemically modulated DUSP6 phosphatase activity. 
Reactive oxygen species (ROS), produced as a result of 
treatment of cells with agents such as H 2 0 2 , block MKP 
enzyme activity (55,56), thereby resulting in high levels of 
Erkl/2 phosphorylation (Figure 5C, left panel). T47D-YB 
cells treated with either 1 mM H 2 0 2 or vehicle alone, 
followed by R5020, exhibited similar levels of PR-B 
Ser81 phosphorylation (Figure 5C, right panel) despite 
effective DUSP6 enzyme (i.e. phosphatase) inhibition, as 
measured by increased phospho-Erkl/2. DUSP6 protein 
levels remained unchanged in the presence of high ROS 
(Figure 5C, right panel). Phosphorylation on other 
selected PR sites (Ser345 and Ser400) was relatively in- 
sensitive to H 2 0 2 treatment (Supplementary Figure S2). 
These data suggest that DUSP6 enzyme activity is not 
required for PR-B Ser81 phosphorylation, as phospho- 
Ser81 levels remained unchanged even under conditions 
where DUSP6 phosphatase activity was greatly dimin- 
ished. Cumulatively, these data suggest that the DUSP6 
protein, but not its phosphatase activity, is required for 
efficient PR-B Ser81 phosphorylation, indicating that 
DUSP6 serves as a scaffolding protein that supports 
ck2-dependent PR-B Ser81 phosphorylation. 

PR-B Ser81 phosphorylation is required for STAT5A 
and Wntl expression 

To link CD domain-dependent regulation of PR-B Ser81 
to gene expression, we examined known PR isoform- 
specific target genes for sensitivity to disruption of the 
CD domain. STAT5A (signal transducer and activator 
of transcription 5A) (25,57) and Wntl (wingless-type 
MMTV integration site family, member 1) (41,58-60) 
are primarily regulated by PR-B in response to progestin 
(25,41). To study the effects of PR-B Ser81 on STAT5A 
and Wntl gene expression, we used a previously well- 
characterized PR-B phospho-mutant that cannot be 
phosphorylated on Ser81, S79/81A PR-B (Ser81 and 
nearby Ser79 were mutated to ensure that this PR 
mutant remained fully resistant to phosphorylation on 
Ser81) (19). T47D cells expressing empty vector control 
(PR-null) or wt, mCD or S79/81A PR-B were treated 
with progestin (0, 1 and 6h), and mRNA was harvested 
for RT-qPCR analyses. After progestin treatment, cells 
stably expressing wt PR-B robustly activated STAT5A 
and Wntl relative to PR-null cells, whereas cells express- 
ing S79/81A PR-B exhibited greatly diminished STAT5A 
and Wntl mRNA relative to cells expressing wt PR-B 
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Figure 5. PR-B's interaction with DUSP6 is required for PR-B Ser81 phosphorylation. (A) DUSP6 is necessary for ligand-induced PR-B Ser81 
phosphorylation. Left: T47D cells stably expressing wt PR-B (T47D-YB) were transfected with 50 nM nonsilencing (NS) or DUSP6 siRNA. Seventy- 
two hours after transfection, cells were treated with lOnM R5020 or EtOH for 60min. Lysates were analyzed via western blotting using p-S81. PR-B 
and DUSP6 antibodies. Right: Densitometry was used to determine the ratio of Ser81 phosphorylated PR-B in R5020-treated cells to total PR-B in 
cells transfected with NS or DUSP6 siRNA. (B) Knockdown of DUSP6 increased MAPK activation. T47D-YB cells were transfected and treated as 
in (A), and lysates were analyzed via western blotting using p-Erkl/2 and Erkl/2 antibodies. (C) DUSP6 phosphatase activity is not needed to 
mediate ligand-induced PR-B Ser81 phosphorylation. Left: Schematic showing how ROS production (increased through treatment with H 2 0 2 ) 
decreases DUSP6 phosphatase activity, which subsequently leads to an increase in Erkl/2 activity and phosphorylation. Right: T47D-YB cells 
were pretreated with 1 mM LL0 2 for 20min, followed by lOnM R5020 (or EtOH) for 30min. Lysates were analyzed via western blotting using 
P-S81, PR, p-Erkl/2, Erkl/2 and DUSP6 antibodies. 



(Figure 6A). Interestingly, cells expressing mCD PR-B dis- 
played an intermediate (reduced expression) phenotype, 
consistent with our finding that this mutant is only 
weakly phosphorylated on Ser81 over time (Figure 3B). 
STAT5B expression remained unaffected in cells express- 
ing wt or mutant receptors (Supplementary Figure S3A). 
Notably, cells expressing the S79/81A phospho-mutant 
PR-B were phenotypically identical to cells expressing wt 
PR-A with regard to STAT5A and Wntl gene expression 
(Supplementary Figure S3B), confirming that these are 
PR-B-specific target genes. Additionally, tissue factor 
(TF) mRNA expression, which is independent of 
PR-B Ser81 (19), was similar in cells expressing either 
wt or S79/81A PR-B (Supplementary Figure S4A). 
Taken together, these data suggest that PR-B Ser81 phos- 
phorylation — which is dependent on a CD domain- 
mediated scaffolding interaction between PR-B, DUSP6 



and ck2 — is required for expression of select PR-B target 
genes (STAT5A and Wntl) known to be critical mediators 
of mammary gland development, stem cell self-renewal 
and breast cancer cell proliferation. 

PR-B CD domain is required for JAK/STAT-dependent 
transcriptional responses 

GSEA is a powerful computational tool that can be used 
to determine whether publicly available gene sets (via the 
MSigDB) are significantly enriched in gene expression 
data sets (39,40). GSEA can identify gene sets that are 
significantly regulated in a particular microarray sample 
group (e.g. wt PR-B or mCD PR-B). Using GSEA, we 
compared ligand-regulated wt and mCD PR-B expression 
data sets (wt PR-B R5020/EtOH versus mCD PR-B 
R5020/EtOH) and identified enriched gene sets from the 
c5 Gene Ontology collection. Interestingly, genes from the 
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Figure 6. Endogenous phospho-Ser81 PR-B-target genes are JAK/STAT-dependent. (A) Select genes regulated by PR-B Ser81 phosphorylation. 
T47D-Y cells stably expressing either wt PR-B, mCD PR-B, S79/81A PR-B or unmodified (PR-null) cells were starved for 18 h in serum-free media, 
followed by treatment with 10 nM R5020 or EtOH for 0-6 h. mRNA levels were analyzed by qPCR. Error bars represent ±SD. (B) GSEA 
comparing gene expression data sets obtained using the Illumina Microarray described in Figure 3. Left: wt PR-B R5020/EtOH compared with 
mCD PR-B R5020/EtOH is shown as compared with the c5 Gene Ontology database, JAK_STAT_Cascade gene set. Right: Leading Edge analysis 
within GSEA (c2 Curated Gene Set) included 229 gene sets significantly enriched (FDR < 0.05) in wt PR-B cells, as compared with mCD PR-B. A 
subset of 38 gene sets (boxed region) contained substantial overlap (green shading) and were significantly involved in interferon signaling pathways. 
See Supplementary Table SI for gene set names and statistics. (C) Ser81 -dependent PR-B target genes are JAK/STAT-dependent. T47D-Y cells 
stably expressing wt PR-B (T47D-YB) were starved for 18 h in serum-free media. Cells were then pretreated (60min) with 50 uM AG490, followed by 
lOnM R5020 for 6h. mRNA levels were analyzed by qPCR. Error bars represent ±SD. (D) PR-B Ser81 phosphorylation is unaffected by JAK/ 
STAT-inhibition. T47D-YB cells were starved for 18 h in serum-free media. Cells were then pretreated (60min) with 50 uM AG490, followed by 
lOnM R5020 for 30min. Lysates were analyzed via western blotting using p-S81 and total PR antibodies. 
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JAK/STAT signaling pathway (JAK_STAT_Cascade 
gene set) were significantly enriched in cells expressing 
wt PR-B relative to cells expressing mCD PR-B 
(Figure 6B, left panel), suggesting that mCD PR-B loses 
the ability to regulate genes in the JAK/STAT pathway. 
Leading edge analysis is a deeper examination used to 
compare only the significant gene sets (identified through 
GSEA) to each other to identify the following: (i) the in- 
dividual genes that are highly associated in a particular 
sample group and (ii) the core gene sets that contain the 
majority of those highly associated genes. These methods 
help identify the pathways (and the core subset of genes 
within these pathways) that are significantly associated 
with a sample group. We performed leading edge 
analysis on gene sets that were significantly enriched 
(229 c2 Curated gene sets with FDR < 0.05) in cells ex- 
pressing wt PR-B relative to mCD PR-B. Our analyses 
identified 38 gene sets that contained substantial overlap 
(Figure 6B, right panel; Supplementary Table SI) and 
revealed a substantial overlap of gene sets involved in 
interferon signaling, which is primarily mediated via the 
actions of STAT proteins. Cumulatively, these data 
suggest that the CD domain in PR-B is critical for 
progestin-dependent regulation of interferon/JAK/ 
STAT-related signaling pathways. 

Interestingly, female STAT5 and PR-B knockout mice 
exhibit similar developmental blocks in mammary gland 
alveologenesis (61). Furthermore, numerous studies have 
implicated the JAK/STAT pathway in PR target gene 
regulation (25,57). Notably, progestin-induced expression 
of selected PR-B target genes, including HSD11P2 and 
STAT5A, was significantly blocked by the JAK/STAT in- 
hibitor, AG490 (62). Similar to STAT5A (Figure 6A), we 
previously showed that progestin-induced expression of 
HSDllp2 requires phospho-Ser81 PR-B (19). To test the 
requirement for JAK/STAT signaling in PR-B Ser81- 
dependent regulation of Wntl expression, T47D cells 
stably expressing wt PR-B were pretreated (1 h) with 
AG490 (50 uM) followed by progestin (10 nM R5020 for 
6h) or respective vehicle controls. AG490 significantly 
diminished progestin-induced Wntl mRNA expression 
(Figure 6C); HSDllf32 and STAT5A were included as 
additional Ser81 -dependent PR-B target genes and 
positive controls (62). In contrast, c-myc, a well- 
characterized PR target gene whose regulation is inde- 
pendent of both PR-B Ser81 (19) and JAK/STAT (62) 
remained unaffected by AG490 inhibition of JAK/STAT 
signaling (Figure 6C). Importantly, progestin-induced PR- 
B Ser81 phosphorylation was unaffected by pretreatment 
of cells with AG490 (Figure 6D). Thus, in addition to PR- 
B Ser81 phosphorylation, JAK/STAT signaling is also 
required for progestin-induced expression of a subset of 
PR-B target genes (STAT5A, Wntl, HSD11(32). 

To determine the role of STAT5 in coregulating 
phospho-Ser8 1 PR-B-dependent transcription, we analyzed 
a publicly available (http://cistrome.org/Cistrome/ 
Cistrome_Project.html) PR chromatin immunoprecip- 
itation ChlP-chip data set for the presence of STAT5 
binding sites within or nearby PR binding sites. These 
ChlP-chip data were created using T47D breast cancer 
cells treated with vehicle or estrogen (to increase 



endogenous PR expression) followed by progesterone 
treatment for 45min. Interestingly, CEAS [cis-regulatory 
element annotation system; (63)] analysis revealed a 1.8- 
fold enrichment of STAT5 DNA sequence binding motifs 
within PR binding sites as compared with random genome 
sampling (P = 8.22 x 10 -11 ). These data suggest that 
STAT5 may function as a pioneer factor by 'opening' 
sites in chromatin for subsequent transcriptional activation 
by PR-B, DUSP6 and ck2. 

PR-B Ser81 phosphorylation is required for binding to an 
enhancer region upstream of the Wntl promoter 

Wntl is a potent breast oncogene and stem cell regulatory 
factor/morphogen. We previously demonstrated the exist- 
ence of a PR-driven autocrine loop in which frizzled recep- 
tors, activated by progestin-induced Wntl, transactivated 
epidermal growth factor receptor (EGFR), leading to 
increased cyclin Dl expression and breast cancer cell pro- 
liferation (41). Importantly, knockdown of Wntl blocked 
progestin-induced breast cancer cell growth in soft agar. To 
further understand how phospho-Ser81 PR-B regulates 
Wntl expression in collaboration with STATs, we per- 
formed an in silico analysis of Wntl promoter and 
enhancer regions. We identified four putative full-length 
PRE binding regions (Figure 7A; triangles), including a 
site located in the proximal enhancer region (~ 1 kb 
upstream from the transcriptional start site [TSS]; PRE1) 
and three sites located downstream of the TSS (+44, 65 and 
66 kb from the TSS; PREs 2-4, respectively). To determine 
whether PR-B directly regulates Wntl, we performed ChIP 
assays to detect relative PR-B recruitment to sites within the 
Wntl enhancer region. After cross-linking and sonication 
was performed, lysates from vehicle- (EtOH) or R5020- 
treated T47D-YB cells were subjected to ChIP using PR- 
specific antibodies. PR-null (T47D-Y) cells served as a 
negative control. In the presence of ligand, we detected 
robust recruitment (~ 16-fold) of wt PR-B to PRE1, while 
moderate levels of PR-B recruitment (~ 2-4-fold with 
ligand) were detected at the other PREs located down- 
stream of the Wntl TSS (PREs 2, 3 and 4; Figure 7B). 
We next tested whether DUSP6 and ck2 were present in 
the transcription complexes at PRE1, the same region 
where wt PR-B was detected. Both DUSP6 and ck2 
proteins were IP'd from similarly prepared lysates of 
ligand-treated T47D-YB cells; PR-B was also IP'd as a 
positive control (Figure 7C). In addition, STAT5-specific 
antibodies revealed that STAT5 was present at this 
Wntl regulatory site, although protein recruitment 
was not ligand-dependent (Figure 7C, bottom right). 
Cumulatively, these data demonstrate direct regulation of 
the Wntl enhancer at region PRE1 by PR-B, DUSP6 and 
ck2. Although these molecules were recruited on progestin 
treatment, STAT5 was pre-associated with this site. 

PR-B Ser81 phosphorylation mediates recruitment of 
PR-B to Wntl and STAT5A enhancer regions 

Finally, we examined PR-B recruitment to Wnt PRE1 in 
cells expressing either wt or S79/81A PR-B. In contrast to 
efficient recruitment of wt PR-B to the Wntl enhancer 
region, S79/81A PR-B failed to be recruited to PRE1 at 



8936 Nucleic Acids Research, 2013, Vol. 41, No. 19 



Wntl 



7$ 



A 



PRE 



7T7E~ 

+44 +65-66 



0 PR-null 
□ wt PR 



Distance from TSS (kb) -1 
B 



15 



20 
18 
16 
14 
12 
10 
S 
6 
4 
2 
0 



] 



PRE1 



PRE2 



PRE3 



n 








5 1 K*? 


p*7 







PRE4 



□ EtOH 
■ R5020 




□ EtOH 
■ R5020 



Antibody 




Antibody 



8.0 

< 7.0 

Q 6.0 

t£ 5.0 

Q. 4.0 
C 

— 3.0 

9: 2.0 

~™ 1.0 

0.0 



u 



□ EtOH 
■ R5020 



Antibody: DUSP6 



IgG 



PR recruitment to Wntl PRE1 



< 

Q 
+^ 

=3 
Q. 
c 



20.0 
18.0 
16.0 
14.0 
12.0 
10.0 



PR-B variant: 




□ EtOH 
■ R5020 



S79/81A IgG 



6.0 

< 

Z 5.0 

Q 

4-1 4.0 

2.0 

5 - 




□ EtOH 
■ R5020 



I I 



Antibody: STAT5 



IgG 



PR recruitment to STAT5A enhancer 



< 

Q 



Q. 

C 0.6 




□ EtOH 
■ R5020 



0.0 j 

PR-B variant: Wt 



S79/81A 



IgG 



Figure 7. PR-B regulates Wntl transcription through binding to Wntl enhancer regions. (A) Schematic of the Wntl promoter/enhancer regions. 
PRE sites (triangles) are located upstream and downstream of the Wntl transcription start site (TSS; denoted with an arrow). Distance of the 
regulatory regions from the TSS is listed in kilobases (kb). (B) PR-B is recruited to Wntl PREs. T47D-Y cells stably expressing wt PR-B or 
unmodified cells (PR-null) were serum-starved for 18 h. Cells were then treated with 10 nM R5020 or EtOH for 60min. Fixed lysates were subjected 
to ChIP with antibodies against PR-B or species-specific IgG (control), and qPCR was performed on the isolated DNA using primers designed to 
amplify the respective regulatory region (PREs 1-4). Recruitment of PR-B to PREs is shown as fold change between R5020- or EtOH-treated cells 
(R5020/EtOH). Error bars represent ±SD of triplicate experiments. (C) PR-B Ser81 phosphorylation-scaffolding complex is recruited to Wntl 
enhancer. T47D-Y cells stably expressing wt PR-B were serum-starved for 18 h. Cells were then treated with lOnM R5020 or EtOH for 60min. Fixed 
lysates were subjected to ChIP with antibodies against PR-B, ck2, DUSP6, STAT5 or species-specific IgG (controls), and qPCR was performed on 
the isolated DNA using primers designed to amplify Wntl PRE1. ChIP experiments were performed in triplicate. A representative experiment is 
shown here; error bars represent ±SD of technical replicates. (D) PR-B Ser81 phosphorylation is necessary for PR-B recruitment to Ser81 -regulated 
genes. T47D-Y cells stably expressing wt PR-B or S79/81A PR-B were serum-starved for 18 h. Cells were then treated with lOnM R5020 or EtOH 
for 60min. Fixed lysates were subjected to ChIP with antibodies against PR or species-specific IgG (controls), and qPCR was performed on the 
isolated DNA using primers designed to amplify Wntl PRE1 or a STAT5A-enhancer site. ChIP experiments were performed in triplicate. 
A representative experiment is shown here; error bars represent ± SD of technical replicates. 
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the Wntl enhancer (Figure 7D, left panel). Similarly, 
ChIP assays revealed that wt but not S79/81A PR-B was 
recruited to a regulatory region within the STAT5A gene 
(Figure 7D, right panel). Importantly, both wt and S79/ 
81 A PR-B showed similar levels of recruitment to PRE 
regions associated with expression of TF [Supplementary 
Figure S4B and (19)], indicating that the lack of S79/81A 
PR-B recruitment to Wntl and STAT5A enhancer regions 
is not due to an inherent nuclear localization or DNA- 
binding defect of phospho-mutant S79/81A PR-B. These 
data suggest a possible feed-forward mechanism involving 
phospho-Ser81 PR and STAT5. Cumulatively, these 
results provide a basis for understanding mechanisms of 
PR-B isoform-specific target gene selection and suggest 
that phosphorylation of Ser81 is required for PR-B re- 
cruitment to selected PREs located within JAK/STAT- 
regulated genes. 

DISCUSSION 

Our studies identify a novel protein interaction domain 
(CD domain) in the PR-B N-terminus that mediates inter- 
action between PR-B, DUSP6 and ck2, and is required 
for progestin-induced S-phase cell cycle entry of breast 
cancer cells. DUSP6 binding to PR-B is required for 
ck2-dependent PR-B Ser81 phosphorylation and subse- 
quent regulation of Wntl and STAT5A, two PR-B- 
specific target genes known to be critically involved in 
mammary stem cell maintenance and breast cancer cell 
proliferation. DUSP6 appears to play a novel scaffolding 
role in this complex, as DUSP6 knockdown but not 
inhibition of phosphatase activity blocked PR-B Ser81 
phosphorylation. Cumulatively, these data support a 
model where DUSP6 binds the PR-B CD domain and 
scaffolds activated ck2 to mediate PR-B Ser81 phosphor- 
ylation (54). In addition, we defined Wntl and STAT5A 
as members of the same gene set regulated by phospho- 
Ser81 PR-B and JAK/STAT signaling. This finding 
complements previous studies showing that PR-B Ser81 
phosphorylation is required for progestin-induced expres- 
sion of HSD11(32 and BIRC3 (19) and that both genes are 
sensitive to JAK/STAT pathway inhibition (62). Notably, 
we found that PR-B, STAT5, DUSP6 and ck2 are re- 
cruited to the same hormone-responsive PRE-containing 
region of the Wntl enhancer. Taken together, our data 
demonstrate a role for PR-B Ser81 phosphorylation in 
target gene selection of JAK/STAT-dependent genes that 
are seminal to mammary gland biology (Figure 8). 
Furthermore, our data demonstrate a functional linkage 
between PR-B and STAT5 that provides a paradigm 
for coordinate activation of proliferative gene programs 
in mammary gland development and, pathologically, in 
ER+/PR+ breast cancer. 

CD domains and nuclear receptors 

The CD domain of PR-B is located in the N-terminal BUS 
region of full-length PR-B, a region that is absent from 
other PR isoforms (21). An in silico analysis of the most 
closely related, steroid hormone-activated nuclear recep- 
tors (i.e. ER, glucocorticoid receptor and androgen 
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Figure 8. Model of phospho-PR-B-specific action. PR-B CD domain- 
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genes, such as Wntl. JAK/STAT-dependent regulation of phospho- 
Ser81 PR-B target genes regulates critical genes involved in mammary 
gland development, mammary stem cell maintenance/expansion and 
early events in breast cancer progression. 



receptor) revealed the presence of weak potential candi- 
date CD domain-like regions. However, none contained 
the conserved Erk2 amino acid sequence similarity found 
in the PR-B CD domain. A similar domain was identified 
in PPARy, a more distant member of the nuclear receptor 
family, which mediates an interaction between PPARy 
and MEK1 (64). Interestingly, deletion of the CD-like 
domain in PPARy did not completely abrogate MEK1 
interaction, suggesting that other areas of PPARy also 
make contact with MEK1 or proteins that bind MEK1 
(64). Similarly, in our study, mutation of PR-B's CD 
domain did not completely disrupt DUSP6 binding. 
Additionally, mutation of PR-B's CD domain did not 
alter MEK1 binding (data now shown), indicating that 
the CD domain in PR-B is not required for PR-MEK1 
complex formation (21). Stable PR binding to components 
of MAPK pathways (MEKs or DUSPs) occurs via 
multiple domains (CD domain, p-Pro, ERIDs) in PR 
that are required for robust progestin-induced MAPK 
activity. Moreover, both c-Src and ER are also found in 
protein complexes with components of the MAPK 
pathway (22,23,65). PR interaction with multiple signaling 
molecules illustrates that these pathways are fully 
integrated, and this integrated circuitry provides a basis 
for the highly selective context-dependent regulation of 
PR target genes. Notably, both the CD domain and p- 
Pro region occur in human but not mouse PR-B. This 
fact may explain why mouse models have not strongly 
implicated PR-B as a breast oncogene; PR-B is not well 
expressed in virgin mice relative to PR-A, and mouse PR- 
B lacks key domains required for linkage of mitogenic 
signaling pathways to PR-B-specific gene expression. 

Regulation of ck2-dependent PR-B Ser81 phosphorylation 

ck2 is a ubiquitously expressed kinase with >300 sub- 
strates (16). Unlike traditional growth factor-activated 
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protein kinases that require upstream inputs for full acti- 
vation, ck2 is constitutively active. Although ck2 regula- 
tion is poorly understood, evidence suggests that ck2 
activation is modulated primarily through subcellular lo- 
calization, substrate distribution/complex formation, ck2 
holoenzyme formation, small molecule interactions or 
autophosphorylation (66). In a previous report, we show 
that PR-B Ser81 is phosphorylated by ck2 in the presence 
of progestin, but in the absence of progestin PR-B Ser81 is 
mainly phosphorylated by ck2 during the S-phase of the 
cell cycle, when ck2 is nuclear and exposed to PR-B (19). 
Herein, we discovered that DUSP6 binding via the PR-B 
CD domain provides a mechanism for robust ck2-depend- 
ent phosphorylation of PR-B Ser81. Perhaps PR-B must 
be bound to DUSP6 to 'accept' Ser81 phosphorylation by 
ck2, either because of proximity restrictions whereby 
DUSP6 recruits ck2 into close proximity with PR-B 
Ser81 or because of substrate conformation changes in 
which DUSP6 binding to PR-B induces conformational 
changes that permit ck2-dependent Ser81 phosphoryl- 
ation. Importantly, we observed constitutive PR-B 
binding with DUSP6 in the absence or presence of pro- 
gestins. However, transcriptional complexes containing 
PR-B, DUSP6 and ck2 are clearly recruited to the Wntl 
enhancer (PRE1) in a progestin-dependent manner. 
Whether these protein-protein and protein-DNA inter- 
actions are regulated by additional factors (e.g. growth 
factors or STAT5) or circumstances (e.g. cell cycle) are 
questions requiring further study. Notably, ck2 is 
upregulated in many human cancers, including breast 
cancer (16,20). Preliminary data obtained from a small 
subset of PR-positive breast tumors [described in (13)] 
demonstrated that roughly half contained phospho-Ser81 
PR-B (A. Daniel and C. Lange, unpublished observa- 
tions). These findings suggest that PR-B Ser81 phosphor- 
ylation is clinically relevant, and underscore the 
importance of further study of PR-B phosphorylation 
and associated isoform-specific target gene expression in 
human breast tumors. 

DUSP6 may function as a scaffolding protein to promote 
cancer growth 

DUSP6 is a potent phosphatase responsible for reversing 
Erkl/2 phosphorylation and thus is a negative regulator of 
MAPK activity. Because of its role as a negative regulator 
of MAPK signaling, the central dogma has been that 
DUSP6 functions as a tumor suppressor in cancer, even 
though DUSP6 overexpression was often predictive of 
poor clinical outcomes [reviewed in (67,68)]. Recent 
data, however, have implicated DUSP6 overexpression 
as a negative prognostic marker in cancer development, 
progression and survival in many different types of 
primary cancers and cancer cell lines, including thyroid 
(69), lung (70), myeloma (71), melanoma (72,73), breast 
(74,75), colon (76), cervical (77), pancreatic (78) and glio- 
blastoma (79). These paradoxical data suggest that 
DUSP6 may have multiple biological functions, independ- 
ent of its long-studied role in attenuating MAPK activa- 
tion. Herein, our data support a novel mechanism 
in which DUSP6 functions as a scaffold for assembly of 



transcriptional coactivators that drive tumor growth. We 
found that DUSP6 phosphatase activity was not required 
for PR-B Ser81 phosphorylation, but rather, DUSP6 
acted as a scaffold protein to bridge the interaction 
between PR-B and ck2, thereby bringing ck2 into close 
proximity with its substrate (i.e. PR-B Ser81) (19,54). 
This scaffolding action represents a unique role for 
DUSP6. Indeed, DUSP6's scaffolding function, rather 
than its phosphatase activity, may provide a molecular 
explanation for the growing body of data linking 
DUSP6 overexpression to poor clinical outcomes in 
many different types of cancer, including breast (74,75). 
In addition, classically defined roles of kinases and phos- 
phatases clearly have broader scopes of action, such as 
bridging pathways previously thought to be unrelated 
(i.e. PR-B regulation by ck2) and direct participation 
in gene regulation as part of transcription complexes 
(19,21,80). 

Functional linkage of STAT5 and PR-B signaling 

Whole genome (cistrome) analyses have identified 'pioneer 
factors' for nuclear receptors. Pioneer factors are 
specialized subsets of transcriptional coregulators that 
bind to transcriptional enhancers, making them compe- 
tent for subsequent transcriptional activation by transcrip- 
tion factors, such as steroid hormone receptors [reviewed 
in (81,82)]. For example, multiple pioneer factors 
(FOXA1, PBX-1, GREB1, AP2-y) have been identified 
for modulation of ER binding, and similar factors have 
been identified for other nuclear receptors (81,83-86). 
Importantly, pioneer factors bind DNA before activation 
of transcription and function to 'open' sites in chromatin 
for subsequent transcriptional activation. We found that 
JAK/STAT pathway inhibition blocked the expression of 
multiple phospho-Ser81 PR-B target genes. We also 
identified STAT5 DNA sequence binding motifs within 
PR binding sites and STAT5 protein associated with 
PR-B target genes. Thus, it is tempting to speculate that 
STAT5 may act as a pioneer factor for phospho-Ser81 
PR-B binding by 'opening' the enhancer regions of 
phospho-PR-B target genes. Interestingly, we also 
observed that phospho-Ser81 PR-B was required for 
STAT5A mRNA expression, suggestive of a feed- 
forward regulatory loop wherein many PR-B and 
STAT5 genes are coordinately regulated. 

Recent reports provide some insight into how PR-B/ 
DUSP6/ck2-containing protein complexes and STAT5A 
coordinate gene expression during breast cancer develop- 
ment and early breast cancer progression. STAT5A and 
Wntl have recently been implicated in PR control of 
mammary stem cell maintenance and mammary gland 
biology (60,87). Like PR-B, STAT5A is required for 
mammary gland development, and both STAT5A and 
PR-B knockout mice have similar defects in mammary 
gland development (87). Progesterone is a known activator 
of STAT5A mRNA and protein expression (25,57); 
however, the mechanism by which progesterone induces 
STAT5A expression is not well understood. Similarly, 
Wnts are important mediators of progesterone action in 
the normal (58) and pregnant mammary gland (59). 
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Previous work published from our lab showed that 
phospho-PR-B-dependent upregulation of Wntl is 
required for breast cancer cell soft-agar growth in 
response to progestins (41). Wnts have recently been 
shown to be critical paracrine mediators of progesterone- 
induced expansion of mammary stem cells (60), and Wnt 
loss as a result of early parity has been linked to protection 
from breast cancer (88). Deregulation of the Wnt/P-catenin 
signaling pathway has been found in many human cancers, 
including breast cancer (89). Interestingly, unlike most 
other cancers, direct mutations of positive and negative 
regulators of the Wnt/P-catenin signaling pathway are 
rarely seen in breast cancer, despite the clear upregulation 
of downstream pathway endpoints, such as (3-catenin 
stabilization and nuclear accumulation (90,91). Notably, 
progesterone/PR-B is a direct activator of this pathway 
(58,60). Potential involvement of these key mediators 
(Wnts, STAT5A) of mammary gland biology in proges- 
tin-induced breast cancer development or early tumor pro- 
gression underscores the need to understand precisely how 
PR-B regulates these genes. 

PR-B Ser81 phosphorylation is a major determinant of 
PR isoform specificity 

We showed previously that phosphorylation of Ser81 is a 
significant determinant of PR isoform specificity; 
mutation of this residue in PR-B (or mutation of CD) 
confers PR-A-specific behavior with regard to target 
gene expression (19) and cell cycle entry (Figure 2C). 
Coordinate regulation of PR-B-specific target genes by 
phospho-Ser81 PR-B and STAT5 may explain the require- 
ment for both factors during the same stage of mammary 
gland development. Notably, paracrine factors (i.e. Wnts, 
RANKL) derived from PR-B-positive progenitors or 
luminal precursor cells are believed to induce self- 
renewal of PR-null stem cells (60,92). PR-A and PR-B 
are usually coexpressed in the same tissues; cells that 
express only a single PR isoform are rare (93-95). A 1:1 
ratio of PR-A to PR-B seen in normal tissues is often 
altered in malignant breast tissues, suggesting that 
balanced isoform action is crucial to normal adult 
mammary gland biology (93,96). PR-A, but not PR-B, 
gene silencing via promoter methylation was significantly 
associated with tamoxifen-resistant breast cancer (97). 
Understanding the critical differences between PR-A and 
PR-B-dependent gene regulation as linked to DUSP6-de- 
pendent PR-B Ser81 phosphorylation by ck2 and JAK/ 
STAT signaling may allow for highly selective isoform- 
specific therapies. Restoration of the balance between 
PR isoform actions may provide an innovative and com- 
plementary approach to existing endocrine therapies. 
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